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Abstract 

We derive a new class of exact solutions characterized by the Szekeres- 
Szafron metrics (of class I), admitting in general no isometries. The source 
is a fluid with viscosity but zero heat flux (adiabatic but irreversible evo- 
lution) whose equilibrium state variables satisfy the equations of state of: 
(a) ultra-relativistic ideal gas, (b) non-relativistic ideal gas, (c) a mixture 
of (a) and (b). Einstein's field equations reduce to a quadrature that is 
integrable in terms of elementary functions (cases (a) and (c)) and elliptic 
integrals (case (b)). Necessary and sufficient conditions are provided for 
the viscous dissipative stress and equilibrium variables to be consistent 
with the theoretical framework of Extended Irreversible Thermodynam- 
ics and Kinetic Theory of the Maxwell-Boltzmann and radiative gases. 
Energy and regularity conditions are discussed. We prove that a smooth 
matching can be performed along a spherical boundary with a FLRW 
cosmology or with a Vaidya exterior solution. Possible applications are 
briefly outlined. 
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1 Introduction 



The relativistic ideal gas' 1 ' is the main conceptual basis for the study of 
matter sources in cosmology' 2 ' 3,4 ' , either in early universe models or in 
classical conditions. The various stages of the thermal history of the uni- 
vese can be described as the evolution of mixtures of various types of ideal 
gases, originally very relativistic and coupled in near thermal equilibrium, 
later becoming non-relativistic and decoupling. A gross simplification of 
this process would group cosmic matter in two main components: "radia- 
tion" (an ideal gas of massless particles) and "matter" (a non-relativistic 
ideal gas). Early universe conditions are characterized by "radiation" as 
the dominant source, while the actual universe is thought of as "mat- 
ter" dominated, with "radiation" being present in the form of relics of 
photon, neutrino and other particles' gases at a very low equilibrium tem- 
perature. 

In order to examine this near equilibrium thermal evolution of ideal gas 
mixtures, conventional wisdom among cosmologists has tended so far to 
use mostly FLRW metrics and perturbations of them' 5 ', usually assuming 
a perfect fluid momentum-energy tensor, with state variables either ob- 
tained from equilibrium Kinetic Theory distributions' 2,3 ' 4 ' , or satisfying 
the so-called "gamma law" (a linear relation between equilibrium pres- 
sure and matter-energy density: p — (7 — l)p). From the perspective of a 
gamma law, the transition from radiation to matter dominance is loosely 
ascribed to the change of the proportionality factor 7 from 7 = 4/3 (ra- 
diation) to 7 = 1 (dust). The latter is widely used as a model of present 
day cosmological matter, since the pressure of non-relativistic matter is 
truly negligible. 

However, these conventional descriptions of cosmological matter sources 
have the following important weak points: (1) The gamma law obscures 
the fact that the actual matter model being considered is that of ideal 
gas mixtures. Instead, the gamma law equation of state with constant 
7 is merely a mathematical construct (with the exception of 7 = 4/3). 
(2) This reasoning applies to dust (7 = 1): the real physical model of 
dominant matter content in a late universe is a non-relativistic ideal gas. 
Dust is merely an approximated description of such a gas at very low 
temperatures. (3) The perfect fluid description (either with a gamma 
law or with Kinetic Theory equilibrium distributions) necessarily assumes 
thermal equilibrium (quasi-static reversible processes). Although early 
cosmological gas mixtures are supposed to be in near equilibrium condi- 
tions, the small deviations from equilibrium can be extremely important' 2 ' 
to account for and explain interactions between particles, cosmological nu- 
cleosynthesis and structure formation. Besides this fact, we should bear 
in mind that dissipative transport phenomena are (in general) present in 
ideal gases, even classical ones at room temperatures ' 6 ' . 

Therefore, bearing in mind the limitations of conventional treatment of 
cosmological matter sources, we propose in this paper to drop the gamma 
law, together with perfect fluid and dust constructs, and to return to the 
original physically sound matter models: ideal gases allowing for dissi- 
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pative fluxes. This requires one to find exact solutions with "imperfect" 
fluid sources in which equilibrium variables satisfy ideal gas equations 
of state, and where non-equilibrium variables are consistent with a posi- 
tive entropy production law and causal transport equations of irreversible 
thermodynamics' 1 ' 7,8 ' 9 ' . 

A recent paper' 10 ' illustrates a possible strategy to find exact solutions 
with an ideal gas fluid model as matter source. The idea is simple: (1) 
consider the metric ansatz of an exact solution usually associated with 
a dust source, (2) replace the dust source with an imperfect fluid, (3) 
impose the equation of state of an ideal gas on the equilibrium state vari- 
ables, and finally: (4) verify that non-equilibrium variables are compatible 
with causal thermodynamics of irreversible processes. This procedure was 
succesfully applied to the spherically symmetric Lemaitre-Tolman-Bondi 
(LTB) metric ansatz' 11 ' 12 ' and to the equation of state of a non-relativistic 
monatomic ideal gas. The resulting models presented two important limi- 
tations: (a) shear viscosity is the only dissipative agent, bulk viscosity and 
a heat flux are zero (the latter necessarily vanishes for LTB metrics), and 
(b) the presence of non-zero pressure gradients with zero 4-acceleration. 
The lack of bulk viscosity is not a problem (it is negligible for a either a 
non-relativistic or an ultra-relativistic ideal gas' 7 ' 8,9 '), but the lack of heat 
flux (with non-zero shear viscosity) makes it necessary to assume that the 
fluid model is limited to a gas evolving along adiabatic but irreversible pro- 
cesses. Fortunately, the above mentioned limitations were compensated 
by the obtention of exact inhomogeneous solutions with thermodynamical 
consistency: the shear viscous pressure satisfies the most general transport 
equation' 9,13 ' (when shear viscosity is the only dissipative flux) with phe- 
nomenological coefficients consistent with those of a Maxwell-Boltzmann 
gas in the non-relativistic limit. 

In the present paper we attempt to generalize the above mentioned 
strategy, to a more general class of metrics and to encompass other ideal 
gas equations of state. Among the most general metric ansatz for a non- 
rotating dust source we have the class of metrics associated with the fa- 
mous Szekeres solutions' 11 ' 12,14 ' (of class I and II), Petrov type D metrics 
admitting (in general) no isometries. In particular, the metric of LTB dust 
solutions' 11,12 ' mentioned above is the spherically symmetric subcase of 
Szekeres class I dust solution. Both classes of Szekeres dust solutions were 
generalized by Szafron' 11 ' 12,15 ' to allow for a perfect fluid source character- 
ized by a geodesic (but expanding and shearing) 4-velocity. This perfect 
fluid generalization has been called the "Szekeres- Szafron" (class I and II) 
of solutions' 11 ' . We shall consider the general metric ansatz of the class I 
Szekeres-Szafron solutions, and proceed as mentioned above by replacing 
their usual dust (or perfect fluid) source with a non-perfect fluid where the 
equilibrium state variables satisfy a generic equation of state comprising 
the following cases: (1) an ultra-relativistic ideal gas, (2) a non-relativistic 
monatomic ideal gas and (3) a mixture of (1) and (2). Again, the obtained 
solutions present the limitations mentioned previously (lack of heat flux, 
with shear viscosity being the only dissipative agent). However, the lack 
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of heat flux can be justified for the case (1) and (3) (in the relativistic 
limit) by the fact that, at high temperatures, shear viscosity is the domi- 
nant dissipative agent' 16 '. The models are mathematically simple and the 
conditions for a consistent thermodynamical intepretation are very sim- 
ilar to those already provided for the LTB case with a monatomic ideal 
gas' 10 ' . We believe these models are a convenient and physically consis- 
tent improvement on existing matter models based on dust and/or perfect 
fluid (with gamma law or not). Finally, the models are ideal theoretical 
tools to explore the rich phenomenology associated to the presence of in- 
homogeneous material sources out of thermodynamical equilibrium. 



2 The relativistic ideal gas. 

A one component relativistic ideal gas is characterized by the equation of 
state' 11 

p = mcnT(fl) - nkT, p = nkT (1) 

(P> ~ K 2 {py 1 kT 

where p, n, T, p are matter-energy and particle number densities, absolute 
temperature and equilibrium pressure, m is the particles' mass, k is Boltz- 
mann's constant and K 2 ,Kz are second and third order modified Bessel 
functions of the second kind. The generalization to an N-componcnt mix- 
ture of relativistic ideal gases is straightforward 



E 



PA, pa = mAC 2 riAT(/3A) — nAkTA, (3a = (2a) 

Kl A 



A=l 

N 

p = ^PA, pA = n A kT A (2b) 

A=l 

where an important particular case occurs if there is local thermal equi- 
librium among the components, so that Ta~T for all A. 

Despite the difficulties in obtaining exact solutions of Einstein's equa- 
tions complying with (1), this equation of state simplifies in the two ex- 
tremes of the temperature and energy spectra: the ultra-relativistic (UR) 
and non-relativistic (NR) regimes, characterized by f3 <C 1 and /3 >• 1 
respectively. Considering the behavior of T(/3) for these cases, we have 

1 + !+0(/3 3 ), p^3nkT (3a) 
/3»1, r^l + ^+0(/T 2 ), p^mc 2 n+^nkT (3b) 
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Another interesting limiting case is that of a binary mixture (A — 1,2) 
where one component is ultra-relativistic and the other non-relativistic (a 
"mater- radiation" mixture), with the pressure of the the UR gas much 
larger than that of the NR one, but comparable to the rest mass of the 
latter. This is obtained by assuming 

/3i < 1, pi w 3mfeTi 
2 3 

pi > 1, P2 « m 2 c n 2 + -U2kT 2 

together with 

U'ikT2 <§C nikTi w m2C 2 n2 

so that 

p w m 2 c 2 n 2 + 3nifcTi, p « mfcTi (3c) 

In particular, if there exists local thermal equilibrium between the two 
components (a very frequent assumption), we have Ti = T 2 = T and so 
the applicability of this approximation strongly depends on the ratio of 
particle number densities m, ri2. For photons and non-relativistic baryons 
the ratio 712/ni « 1CP 9 is a small number, hence for the temperature 
range 10 3 K < T < 10 5 K, characteristic of the radiative era, we have 
nikT <C mkT together with mkT « mc 2 n 2 , and so this approximation to 
a matter-radiation mixture is reasonable' 3 ' 16 '. If there is no local thermal 
equilibrium (for example, after decoupling), then Ti 7^ T 2 , but we can still 
justify (3c) as long as nikTi <C niATi holds. 

We shall describe the three cases (3) through a generic equation of 
state characterized by the ansatz 

2 nkT , , , 

p — mc n -\ -, p = nkT (4) 

7-1 

where 7 is a positive constant, coinciding in the ultra and non-relativistic 
limits of (1) with the ratio of heat capacities at constant volume and pres- 
sure, which in these limits takes the constant asymptotic values 4/3, 5/3. 
With this ansatz we can recover (3a-c) through the following particular 

cases 

4 

Ultra-relativistic ideal gas : m = 0, 7 = - (5a) 

o 

4 

Matter-radiation mixture : m > 0, 7 = — (5b) 

5 

Non-relativistic monatomic ideal gas : m > 0, 7 = - (5c) 

where, in the case (5b), it is necessary to bear in mind that n multiplying 
mc 2 is TI2 , while that multiplying kT is m . However, we shall not use these 
subindices, unless the distinction between 711,712 is not clearly understood 
from the context of the discussion. Also, for classical ideal gases, we 
can generalize (5c) to 7 = 1 + 2/q, for q a natural number, describing 
molecules with q degrees of freedom' 17 ' (the case (5c) corresponds to q = 3, 
a monatomic gas). We show in the following sections that Einstein's field 
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equations can be integrated for Szekeres-Szafron metrics of class I whose 
source is a fluid matter tensor complying with the generic equation of 
state (4). 



3 The Szekeres-Szafron metrics. 

The class I of perfect fluid Szekeres-Szafron solutions is characterized by 
the following metric 111 ' 12 ' 14 ' 151 

ds 2 = -c 2 dt 2 + {uhf (Y'Y dz 2 + Y 2 [dx 2 + dy 2 ] (6a) 

Y = j- (66) 

where 3> = §(t,z), h — h(z), Y' = F z and the function v = v(x,y,z) is 
given by 

v = A(z)(x 2 + y 2 ) + B 1 (z)x + B 2 (z)y + C(z) (6c) 

with A,Bi,B2,C being arbitrary functions. This metric admits, in gen- 
eral, no isometries, but contains as particular cases spherically, plane and 
pseudo-spherically symmetric solutions (3-dimensional isometry groups 
acting on 2-dimensional orbits), as well as FLRW spacetimes with posi- 
tive, negative or zero curvature of spacelike slices' 11 ' . 

We shall consider as matter source for (6) a viscous fluid characterized 
by the stress-energy tensor 

Tab a b , iab , ab /~\ 

— pUU+ph+TT (7) 

7 ab ab . —2 a b b f\ a ri 

h — g + C U U , TTabU — 0, 7T a = 

associated with the balance and conservation laws 

T ab . b = (8a) 

(nu a ). a = (86) 

[nsu a ), a > (8c) 

where s is the entropy per particle. 

In the comoving representation (u a = cS a t) for (6)- (7), the 4-acceleration 
vanishes (u a = u a - b u b — 0) and the remaining non-zero kinematic in- 
variants are the expansion scalar: 6 = u a :a and shear tensor: o ab = 
u (a;b) ~ (1/3)6 h a t- These invariants, together with the shear viscous 
pressure (a traceless symmetric tensor like a a (,), are given by 

Y' ?Y 

e =F7 + |- 0) 

a\ = diag {0,-2<t,<t,<j}, a = i - ^ (10) 
7r a b = diag [0, -IP, P, P] (11) 
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where Y = u a Y >a = Y t — $,t/f and P — P(t,x,y,z) is determined by 
the field equations. 

The integration of the conservation equation (8b) for (6) leads to 

(12) 



Y 2 Y> 

where N = N(x, y, z) is an arbitrary function (the conserved particle 
number distribution). For the mixture in (5b), we have assumed the two 
components to have independently conserved particle number densities 
and to be characterized by the same 4-velocity, hence: N = iVi + N2 in 
this case. Einstein field equations for (6) and (7) yield 

(13a) 





[Y(Y 2 + Kc 2 )}' 


YW 


[Y(Y 2 + Kc 2 ) + 2YY 2 




3Y 2 Y' 


Y 


'Y 2 + Kc 2 +2YY~ 


6Y' 


Y 2 



«P=~ — ^T, (13b) 



(13c) 



where k — SttG/c 2 and the function K — K(x, y, z) is given by 

Ah 2 (AC - B\ - Bl) - 1 
K = ^ H 2 » 2 ) ^ 
This function' 11 ' 15 ' is related to the curvature of the hypersurfaces t = 
const, everywhere orthogonal to u a . The energy and momentum balance 
(8a) for (7) is given with the help of (4) by the following equations 

p + 7p9 + 6Po- = (15a) 
6Py + (p + 2P)' = (156) 

/ ,x \ ' ,y 

where v(x, y, z) is given by (6b). We have written equations (13) in terms 
of Y rather than in terms of $ and v, thus obscuring the dependence of 
the metric on the coordinates (x,y). However, equation (15c) places a 
strong restriction on the way in which initial value functions and state 
variables may dependend on these coordinates. This is illustrated in the 
following section. 



4 Integration of the field equations. 

Considering the field equations for (6) and (7) satisfying (4), we look 
first at the possibility of a reversible evolution along all fluid worldlines, 
associated with a perfect fluid tensor (ir a b = in (7)), and then examine 
the irreversible case with shear viscosity. 
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4.1 Perfect fluid case. 



Szekeres-Szafron perfect fluid solutions of class I follow by setting P = 
in (13c), which together with (13b), yields 

Y 2 +Kc + 2YY = -np{t)Y 2 , (16) 

where p(t) is the equilibrium pressure, an arbitrary time-dependent func- 
tion that must be specified in order to have determined solutions. Since p 
obtained from (16) and (13a) depends also on the spatial coordinates, it is 
an obvious (and a known) fact' 11 ' 15 ' that these perfect fluid solutions are 
incompatible with the so-called "barotropic equations of state" of the form 
V = P(p)- However, (4) is not barotropic, but a two-parameter equation of 
state expressible as p = p(n,T) and p — p(n,T), and so it is not obvious 
whether these solutions admit this equation of state or not. It is also 
known' 18,19 ' that the class I of Szekeres-Szafron solutions, whose source is 
a one-component perfect fluid, is compatible with a two-parameter equa- 
tion of state and with the integrability of the equilibrium Gibbs equation 
only if it admits an isometry group acting on orbits of two or more dimen- 
sions. Therefore, possible perfect fluid Szekere-Szafron solutions of class I 
compatible with (4) must necessarily present this type of symmetry. The 
case of a fluid mixture is not subjected to this restriction, but will not be 
discussed. 

Assuming the perfect fluid source (P = 0) of a Szekeres-Szafron metric 
to be compatible with the integrability of the equilibrium Gibbs equation, 
equations (15a) and (16) imply = 0(t), but then (8b), (rewritten as: 
h + n0 = 0), yields n = n(t). From (4) and (12), we have p — p(t) and 
T = T(t), while (9) implies in this case $(f, z) = $i(i)$2(z), all of which 
means that the metric reduces to that of a FLRW spacetime. Therefore, 
we conclude that a Szekeres-Szafron fluid solution, complying with (4) 
and not reducible to a FLRW metric, necessarily requires a momentum- 
energy tensor with non- vanishing anisotropic pressure (or shear viscosity). 
In other words: reversibility along all fluid worldhnes is only possible for 
the FLRW subcase of (6) (see section (9.3)). 



4.2 The viscous fluid case. 

If P =/= 0, the field equations (12), (13a) and (13b), subjected to the 
equation of state (4), still provide a closed system and, consequently, 
Einstein's field equations can be solved independently of (13c). In this 
case, equation (12) plus the field equations (13a) and (13b) subjected to 
the equation of state (4), yield the constraint 

[(3 7 - 2)Y(Y 2 + Kc 2 ) + 2YY 2 ] ' - 3k (7 - 1) mc 2 N = (17a) 
which integrates to 

(3 7 - 2) (Y 2 + Kc 2 ) Y + 2YY 2 - 3k (7 - 1) M = f(t, x, y), (176) 
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where f(t, x, y) and M(x, y,z) = J Ndz are arbitrary functions. However, 
these functions are not totaly unrestricted and, as shown in section 9.2, 
a consistent dust limit associated with p = P = requires / = 0. Under 
this assumption the specific functional dependence of M on (x, y) is given 
by 

M=^fi. (18) 
v 6 

This follows from consistency of (17) with (6b) and (14) (once we have 
set / = 0), and agrees with (15c). Integrating (17b) again, for / = 0, we 
get 

(Y 2 + Kc 2 )Y = KM+ y ^ J (19) 

where J = J(x, y, z) is an arbitrary integration function. In order to infere 
the form of this function, we substitute Y + Kc from (19) into (13a) 
and (using (12) and (18)) compare it with p given by (4), resulting in 

NkT J . . 

j — i y3(7-i) v ' 

Therefore, by analogy with (18), we can define J in terms of an initial 
value for the function at the left hand side of (20a). This leads to 

J = .UoY^'-V, U = = / ^ dz (20b) 

v 3 J 7 - 1 

so that 

where, as with (18), the restriction on the dependence of Uo (and thus, 
To) on (x, y), follows from (6b), (14) and (15c). From now onwards, a 
subindex"0" will denote the "initial value" of the function it is attached to. 
That is, the function evaluated along an arbitrary spacelike hypersurface 
orthogonal to u a and marked by the initial time t = to- Using hqYqYq = 
N from (12), we can rewrite (18), (19) and (20b) as 



Y 2 = — 



M + Uo 



'F n X 3(7-1) 



- Kc 2 (21a) 



® 

\t 2 f vV'^ rr f nokT Y 2 Y ' J 
M = mc I noY Y dz Uo = I dz 



(216) 



in terms of initial value functions no, To, and K, with clear physical and 
geometric meaning. However, no, To are not completely general functions 
since they must satisfy the restrictions given by (18) and (20b), though 
M, Uo are totaly arbitrary and the function v contains five extra arbitrary 
functions of z, so these restrictions basically affect the dependence of no 
and To on the coordinates (a;, y). Equation (21) is a generalized Friedmann 
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equation (or evolution equation) for fluid layers labelled by the comoving 
coordinates. Again, for the case (5b), the function no (like n and N) 
multiplying mc 2 and kT correspond to a separate mixture component. 



5 The state variables. 



Before proceeding to integrate (21), it is useful to derive semi-determined 
forms of the state variables and kinematic invariants by inserting (21) into 
(12) and (13). This yields 



P = 



n [y 



YqY' Yj/Yp 
Y'/Y 

3(7-1) 



mc no + 



nokT 

7 



kTo (Yo 
-1 \Y 



3(7-1) 



Yq/Yq 
Y'/Y 



(22) 
(23) 
(24) 
(25) 
(26) 



6(7 



y 

Y Y 



men 



1 - 



2YY' ; 2 

— — 6 = k [mc no 



Y Y 
3U 



3M Y' 

(¥)' 



2P 



XL 

Y 



Kc 2 
Y 2 



K 



2Y' \ 

Y ) 

(27) 



3M Y' 
+ ~M>Y 



nokT (Yo 

7- 



37 



Uo 



V 

J- n 



i + TFT (7-1)^ + (2 



Y 2 



Y f-l^' : (2S) 



4Y^ 



where 



* = 1+3(7- 1; 



Uo 
U> 



XL 

Yo 



XL 

Y 



!lsl| 



Uo 
U' 



3(7 - 1) 



XL 
Yo 



37 



y 



(29a) 



(296) 



The new solutions, characterized by (21)-(29), become fully determined 
once Y, Y' are found by integrating (21) subjected to initial conditions 
given by the initial value functions no, To, Yo and K. Because of (6b), (6c), 
(18) and (20b), the initial value functions are really the z dependendent 
functions A, Bi, £>2, C, h and M, Uo- Notice that it is always possible to 
eliminate any one of these funcions by rescaling the z coordinate. 
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6 General integral and specific solutions. 



Equation (21) can be transformed into the following adimensional quadra- 
ture 



Y 



(t - to) 



ry 



y(37"2)/2 d y 



where: 



Y_ 

Y 



kUq kUo 



(30) 



kM 



c 2 Y c 2 i/ 2 $o : 



c 2 Yq c 2 i/ 2 <& 



so that (30) can also be though of as the integral yielding $ (which, 
from (6a) yields Y). We integrate (30) explicitly for the three cases of 
interest mentioned in equations (5) . From the results of these integrals it 
is possible to obtain the gradients Y' and to evaluate explicitly (22)-(29). 

6.1 Ultra-relativistic ideal gas. 

Setting 7 = 4/3 and m — (or, equivalently fj, = 0) in (30), we obtain the 
following solutions 



Case K = 



Case K ^ 



Y = Y 



l±^(*-*>) 



1/2 



(31a) 



^-(t-to) =±A" 1 Vu-K-y/u- Ky 2 (316) 

6.2 Matter-radiation mixture. 

This case is similar to the previous one: 7 = 4/3 but /i > in (30) (so, 
M > 0). As mentioned previously, for this case the functions no, n or 
N multiplying mc 2 or kTo correspond to each separate component of the 
mixture. 



Case K = 0. 



§V^^(* - *o) = y/y + 6 {y - 26) - VT+5 (1 - 25) (32a) 

where 

s _ u _ Uo 
~ n~ M 

Case K ^ 0. 



11 



A = In 



Kjiy - Ku - 2Ky + fj, 



-Kh-Klu -2K + [i 



K < 



. u-2K . u — 2Ky 

A = arcsin — ■ arcsm ■ 



K > 



A = /x — 4A^ 



6.3 Non-relativistic monatomic ideal gas. 

This case follows by setting: 7 = 5/3 and /x > in (30). The correspond- 
ing integrals are combinations of elliptic integrals and algebraic functions. 
They have been obtained with the help of standard text books dealing 
with these integrals' 20 ' 



Case K = 0. 



\^Y (* - *o) = Vv (V 2 + S )~ v^+t + ^ (F - F ) (33a) 
where 

_ u U 

and F is the elliptic integral of the first kind with modulus l/y/2 and 
argument ip given by: cos tp = (y — VS)/(y + \/5). (see equation 241.00 of 
reference [20]). 

Case K ^ 0. 

This case follows from equations 260.03 and 341.05 of reference [20]. It 
requires re-writing (30) as 



Y o ' Ji {y(~Ky* + ) iy 2 +cu)} 1/2 



fdy 



± f — m (336) 

h |-^y(y-6)[(y-6 1 ) 2 + a 2 ]} ' 



which integrates to 



«( t _ to)=± _^ 
Y K ' (U + V) 

where 



(t-to)=± — ; rn2 g [A(fre)- A(fo,e)] (33c) 



A(^O=F(0,e)- T7 ^- J 



R (1) (<^)--^R (2) («M 
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R 



(m) , 



r<t> 
Jo (1 



M 

+ Ccostf) m VI - e 2 sintf 



and F(0, e) is the elliptic integral of the first kind with argument <j> and 
modulus e. The latter parameters plus the remaining quantities in (33b) 
and (33c) are 



b = 



a 2 + 4u 2 + 2[m 
6~Ka ' 



bi = - 



a 2 + 4fi 2 



4pa 



YlKoi 



ai = 



Vs(4^ 2 



VIKa 



r2 _ a 4 - 2/ia 3 



U z = 



i 3 a + 16/j 4 



36a 2 



2 = q + 4(i 2 a 2 + 16u 4 
12a 2 



6Ka 


(2/i - a) - ^3(a 2 - 2/ta + 4/i 2 ) 


y- 


(2M- 


a)(a 2 + 2/m + 4/i 2 ) 


6Ka 


(2/i - a) + v/3(a 2 - 2/xa + 4/i 2 ) 


2/- 


(2^- 


a) (a 2 + 2/m + 4/i 2 ) 



9 = 



V + U 




a 3 = 108a; + 8u 3 + 12V3K [lu(27uj + 4u 3 )] 



7 Conditions for thermo dynamical con- 
sistency. 

We have obtained a class of exact solutions whose equilibrium variables 
satisfy the equations of state (3)-(5). The entropy per particle, s, and 
the shear viscous pressure, ir ab , (the only dissipative flux present in the 
matter tensor) were not used in the integration of the field equations. It 
is necessary now to examine the consistency of these quantities within 
the framework of a physically acceptable theory of irreversible thermody- 
namics. Concretely, this means verifying that ir ab and s satisfy a suitable 
transport equation and the entropy balance law (8c), for physically rea- 
sonable phenomenological coefficients compatible with (3)-(5). 

The oldest approach to irreversibility is the so called "classical" theory 
of irreversible thermodynamical processes' 7 ' , based on the hypothesis of 
"Local Equilibrium". This implies considering s as depending only on 
equilibrium variables and obtained from the integration of the equilibrium 
Gibbs equation: Tds = d(p/n) +pd(l/n). For: n, T, p, p related by the 
equation of state (4), this means 



s (e) = s + feln 



1/(7-1) 



no 
n 



= so + k In 



Y'/Y 



1/(7-1) 



(34) 



13 



where so is an initial value function. In order to apply (34) to the case (5b), 
we have assumed a single T for the mixture (local thermal equilibrium 
between the components). In this case, n, no appearing in (34) can be 
considered particle number densities of the ultra-relativistic component. 

It is a well known fact that equation (34), together with (8c), yield 
parabolic transport equations' 7-9 ' , and so the classical theory is applica- 
ble to physical systems whose microscopic time and length scales (mean 
free path and mean collision time) are much smaller than the macro- 
scopic evolution times or characteristic length scales. However, these con- 
ditions are not always applicable to astrophysical and cosmological sys- 
tems, both non-relativistic and relativistic. The former (globular clusters, 
galaxies, dust clouds) are practically collisionless, or have mean collision 
times comparable to their evolution time, the latter (early universe gas 
mixtures) cannot be compatible with transport equations violating causal- 
ity. The alternatives to the classical theory are the theories generically 
grouped under the term "Extended (or causal) Irreversible Thermodynam- 
ics" (EIT)' 7-9 ' . 

EIT assumes that s also depends on the dissipative stresses. For the 
momentum-energy tensor (7), where shear viscosity is the only dissipative 
stress, EIT associates a generalized entropy current satisfying (8c) and 
relating the deviation from equilibrium due to viscosity 

{nsu a )- a > => s > 0, s = s^' — air ab -K ab (35) 

where a is a phenomenological coefficient and s^ e ' is given by (34). If 
shear viscosity is the only dissipative stress, the most general available 
transport equation in the formalism of EIT is' 9 ' 13 ' 



T-kad h c a hf + n ab 



+ 2r t a ab = (36) 



where -k ab = ir ab . c u c , and r, n are the relaxation time and the coefficient of 
shear viscosity, phenomenological quantities whose form depends on the 
properties of the fluid. The transport equation of the classical theory is 
recovered from (36) by assuming r = . Also, within EIT, a simpler or 
"truncated form" without the term involving u c and T in (36) is often 
suggested' 7,8 ', for no other reason than to deal with a mathematically 
simpler transport equation. However, althogh this truncated form corrects 
the problems of acausality and unstability plaging the classical theory, we 
shall consider the full transport equation (36) , since numerical studies with 
FLRW spacetimes' 21 ' and other theoretical arguments' 22 ' reveal adverse 
physical effects, and so, advice caution in using truncated versions of 
viscosity transport equations. 

Equations (35) and (36) couple to Einstein's equations through the 
dependence of state variables on the metric functions. However, the ful- 
filment of these equations cannot be verified until an explicit form for 
the phenomenological quantities a, r, r\ is somehow specified. As a refer- 
ential value of the forms that one should expect for these quantities in 
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dealing with one-component ideal gases associated with equations of state 
like (3)-(5), consider a relativistic Maxwell-Boltzmann distribution near 
equilibrium for which the equation of state (1) holds' 7,8 '. From equations 
(7.1) and (7.10c) of reference [8], and bearing in mind that the factor 
multiplying ir a t is r, we obtain r\ in terms of r as 



V = 



r 2 (/3) 



2/3 2 i + ertpyp 



pr 



(37a) 



where F(/3) is given in equation (1) Regarding the case (5b), the most 
convenient referential system for such a gas mixture is the "radiative 
gas"' 4 ' 7 ', a mixture of photons and non-relativistic particles, for which 
we have: 



r\ = — pr 



(376) 



where p is the pressure of the ultra-relativistic component (photon gas). 
For both cases, (37a) and (37b), the coefficient a is given by 

(37c) 



4t? nT 



Considering the expansion of T(/3) in the ultra-relativistic (/3 <C 1) and 
non-relativistic {fi ^> 1) limits in (37a), we can provide a compact expres- 
sion for T],a comprising these limits (cases (5a) and (5c)), as well as the 
radiative gas (case (5b)) 



7} = b pr, 



1 _ k 
4b pnT 4b P 2 



( 2/3, Ultra- Rel. 
bo = < 4/5, Mixture 
{ 1, Non-Rel. 

(37d) 

where p — nkT was used in (37d). The relaxation time r, for all these 
cases, is given by collision integrals (for example, see eqs (3.41) and (3.45a) 
of [7] for the non-relativistic case). 

The consistency of the solutions with EIT and Relativistic Kinetic 
Theory follows from verifying to what degree the forms of T, p, P and a 
, as given by equations (22)-(29), are compatible with (35), (36) and the 
coefficients (37). Substitution of (21)-(29) into (35), with s^ and a given 
by (34) and (37c) yields 



s = so 



kin 



T \ 1/(7-1) 
To 



3k 
26^ 



so + k In 



Y'/Y 



1/(7-1) 



3k /n 



(38) 



where so(r) is an arbitrary initial value of s. Regarding (36), if we assume 
the Maxwell-Boltzmann form: r\ = bopr given by (37d), with p given by 



1 The quantities (rj, /?2,r,p) here are respectively given by (£s, /32, r\, P) in reference [8] 
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(23), then this transport equation is satisfied for r given by: 

9 U Y'] 2 , fU Y'\ 2 \ 

*-*mY\ +bi {w>r) ) (39) 

While the form of r\ is formally identical to the Maxwell-Boltzmann shear 
viscosity coefficient (37d), and thus needs no justification, the form of r in 
(39) is acceptable as long as this expression behaves as a relaxation param- 
eter for the system, namely: it should be positive for all fluid worldlines, 
and must guarantee that s (computed from (38)) is also positive every- 
where. Ideally, of course, r should also, somehow, approach or relate to 
its corresponding collision integrals from Kinetic Theory of the Maxwell- 
Boltzmann relativistic gas. It must also relate to the mean collision time, 
and so should increse (decrease) as the fluid expands (collapses). Evalu- 
ating s from (37) and comparing with (39), we find that s and r satisfy 



4b () cr 



(nsu a ) 



2rjT 



3k 
bor 



3k /O 
4b T V* 



ny 



(40) 



a relation that also follows directly from (35) and (36). Equation (40) 
implies: s > r > 0. Therefore, since bi > in (39) for the cases 
of interest (6 = 2/3,4/5,1), equations (39) and (22)-(27) lead to the 
following necessary and sufficient conditions for T > 0, r > 0, s > 



* > 



crQ < 



(41a) 
(416) 



while the conditions insuring that s decreases for increasing r (r > <4> 
s < 0) follow from (39)- (40) 



T > 0, 



2a UoY^_ 

*n y 



1 + 3(7-1) 



U^Y 



- T -<o 

T 



(41c) 



The set (41) provides the necessary and sufficient conditions for a theoret- 
ically consistent thermodynamical description of the solutions complying 
with cases (5a-c) of the generic equation of state (4), within the framework 
of equations (34)-(40). The verification of the fulfilment of (41) requires 
the explicit computation of the gradients Y' /Y from (31)-(33). However, 
fulfiling this task in the present paper would make it too long, hence we 
believe it is more convenient to present here the general theoretical out- 
line and leave a proper full examination of thermodynamical features for 
each case (5) in separate papers' 10 ''' 23 '. In particular, it is worthwhile 
to mention that physically reasonable numerical examples satisfying (41) 
have been found' 10 ''' 23 ' for the spherical subcase of (6) and the equations 
of state (5a) and (5b). 
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8 Energy conditions and regularity. 



Together with the conditions for thermodynamical consistency, the solu- 
tions must comply with energy and regularity conditions. From (22)-(29) 
it is evident that Y/Yo =0 marks a scalar curvature singularity charac- 
terized by the blowing up to infinity of all state variables. This singularity 
is analogous to a big bang in FLRW cosmologies, and has been reported 
in Szekeres-Szafron solutions with dust and perfect fluid sources' 24 '"' 25 ' . 
As in these cases, it is not marked (in general) by a constant value of t . 
However, an important requirement for the present models is to select an 
everywhere regular initial hypersurface, so that t = to lies entirely in the 
future (past) of Y/Yo = for expanding (collapsing) configurations. 

Another scalar curvature singularity, of the so called "shell crossing" 
type, occurs if the gradient Y' /Y vanishes. Since this gradient appears in 
the denominator of n, p and P, but does not appear in the expression for 
T, if y'/V = we would have n, p and P diverging but (in general) finite 
T. This is, obviously, a totaly unacceptable situation from a physical point 
of view. Therefore, we must demand the following regularity condition 

y'/y 

YUY Q >° ™ 

to hold along the entire evolution range being considered. Notice that the 
possibility of having an opposite relation sign in (42) is ruled out because 
n would be negative and, for p > 0, it would imply ^ < and so, T would 
be negative. 

The eigenvalues of the momentum-energy tensor (7) are: Ao = — p, 
Ai = p — 2P, A2 = A3 = p + P, hence the weak, dominant and strong 
energy conditions (WEC, DEC and SEC) are given by 



p>0, p + p-2P>0. p + p + P>0, WEC (43a) 

p > 0, \p- 2P\ <p, \p + P\< p, DEC (436) 

DEC plus: p + Sp > 0, SEC (43c) 

With the help of (4) and (22)- (25), and assuming (42) to hold, we can 
translate (43) in terms of the parameters of the solutions, leading to the 
following equivalent set of conditions 



WEC 



mc no 



mc 2 n (J£) 7 + nokT f^y* - ^ j > (446) 



mc 2 no 
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SEC 



DEC 

nokT |* - Q\ < mc 2 n ^ + ^jl* (44d) 

v y o i i — i 

\nokTo !2* + fi| < mc 2 n f ^V' 7 ^ + (44e) 
2 V x o / 7 — 1 

mc 2 n ( z-z- ) + -!—. -n fcT * > (44/) 

V Yo J 7 — 1 

where in the case of the matter-radiation mixture (5b), the initial densities 

no multiplying mc 2 and kT respectively correspond to the non-relativistic 

and ultra-relativistic components. Notice that conditions (41a) and (42) 

are sufficient for the fulfilment of (44a) and (44f), while the fulfilment of 

the remaining conditions can be examined by reducing, with the help of 

(21b), (29a) and (29b), the terms appearing in (44) to the following forms 



3( 7 - 1)U Y'/Y 
n kT Y 3 Y '/Y 



1 3( 7 -l)Pb 
7 - 1 n kT Y 3 



2* + n = 3 



! + 3(7 - 1) 2 W) (, 3 7 -2 F'/y 



2 ^ -* + «- 3 ^ 



7 — 1 7 — 1 



n„fcT„y ( f ^ 3( 7 - 1) y '/y 

" 1+ 3( 7 -l) 2 t/o ^ 3 7 F'/V 



(45c) 



(45d) 



n kT Y a \ 37-iy '/y 0/ 

It is evident that the rhs of (45b) is positive for 7 > 1, while (42) is 
sufficient for the terms in the rhs of (45a) to be positive, and so (41a) 
and (42) are sufficient for the fulfilment of (44b) and (44d). However, 
nothing can be said about (44c) and (44e) until the gradient Y'/Y is 
evaluated. The conditions for thermodynamic consistency do not restrict 
the sign and magnitude of Q, as long as (41b) holds, but conditions (44c) 
and (44d), through (45c) and (45d), do provide further restrictions on 
the gradient of Y controlling how much bigger can |fi| be with respect to 
*, or how large can be \P\/p. Therefore, if \P\/p > 1, then the WEC 
(DEC) might be violated if P < (P > 0). However, this ratio between P 
and p is unphysical since the linear-thermodynamics approach, as that of 
EIT, demands that \P\/p <C 1 because of the near-equilibrium conditions 
on which the theory is based on. Like the fulfilment of conditions (41), 
the fulfilment of (44) for each specific case will be examined in separate 
papers. 



9 Spherical, dust and FLRW subcases. 

The solutions derived and classified in the previous sections contain vari- 
ous important particular cases, obtained either by imposing extra symme- 
tries on (6), by restricting the matter sources or both. For the purpose of 
describing these subcases, it is useful to transform the spatial coordinates 
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of metric (6) to spherical coordinates. This is achieved by the transfor- 
mations: z — r, x = 2tan(#/2) cos(0), y = 2tan(#/2) sin(0), leading to 

ds 2 = -c 2 di 2 + ftV (y') 2 dr 2 + F 2 [d<9 2 + sin 2 8d<f> 2 ] (46) 
Y = ^ = 4A sin 2 (6»/2) + 2 sin 6»(Bi cos + B 2 sin0) + C cos 2 ((9/2) 

where <£> = <£>(£, r), h, A, B\, E>2,C are now functions of r and a prime 
denotes now derivative wrt r. 

9.1 Spherical symmetry. 

From (46), it is clear that spherically symmetric subcases of the solutions 
follow by especializing the functions A, Bi, B2,C to the values 

Si = B 2 = 0, A' = C = 0, 4A = C (47a) 

leading to 

v = \, K = C 2 -^ (476) 

M = M(r), U = U (r) (47c) 

while = /i(r) remains arbitrary. In fact, the deviation of A,Bi,B2,C 
from the values in (47a) "gauges" the deviation of the geometry from 
spherical symmetry. The metric resulting from applying (47) into (46) is 
known as the Lemaitre-Tolman-Bondi' 11 ' 12 ' metric ansatz, and it is usu- 
ally associated with dust solutions. However, this ansatz can also be the 
metric of exact solutions with a perfect fluid source (spherical subcase of 
Szekeres-Szafron class I solutions) or with the type of viscous source exam- 
ined in the present paper. Identifying: $ = Y and h 2 — 1/(1 — fco/ 2 (r)), 
leads to the solution examined in reference [10], where only the equation of 
state (5c) was considered. Other particular cases characterized by isom- 
etry groups acting on 2-dimensional orbits (plane and pseudo-spherical 
symmetries) can be obtained by specializing the functions A, Bi, Bi, C to 
specific constant values (see [11]). 

9.2 Dust subcases and dust limit. 

If no > but T = (or equivalently: M > and Uo =0), we have 
T = p = P = with p = mc 2 n = M'/(Y 2 Y'), leading to the Szekeres class 
I dust solution' 11 ' 12,14 '. If besides these restictions on the state variables 
we include spherical symmetry, we obtain the Lemaitre-Tolman-Bondi 
dust solutions 111,121 . 

However, the solutions must have a consistent dust limit, by which we 
shall mean that setting p — in equations (13) and (15) yields exactly 
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the same results characterizing Szekeres dust solutions or their particular 
cases with isometries. We prove below that, as mentioned in section 4.2, 
a consistent dust limit for the solutions requires f(t,x,y) = in (17b). 
First, we re-write (17b) as 

Y (Y 2 + Kc 2 )+ 2YY 2 = f(t, x, y) + 3( 7 - 1) [kM - Y(Y 2 + Kc 2 )] 
so that substitution into equations (13b) and (13c) yields 

[kM - Y(Y 2 + Kc 2 )] ' 



/tp = -3(7-1) - 



yly, 



P Y ftf 1 \'^^ 1X \ kM-Y(Y 2 + Kc 2 ) ] ,\ 

The condition p = yields 

y i = K,M-g(t,x,y) _ 

which coincides with the Friedmann-like evolution equation of Szekeres 
dust solutions only if g = g(x,y), but then this function can be absorbed 
into M, hence no generality is lost by setting g — 0, leading to 



kP-^ 
6Y 



but from (15a), p — implies P — (or a — leading to FLRW subcases, 
for which P also vanishes), therefore we must have / = for a consistent 
dust limit. 

The solutions are also consistent with the intuitive idea that dust is 
an asymptotic and low temperature limit of an expanding ideal gas. In 
fact, by looking at the Friedmann-like equation (21), it is evident that for 
m > and 7 > 1, this equation tends, for Y/Yq 3> 1, to its equivalent 
evolution equation for dust. Similarly, the form of the state variables 
(22)-(29) reveals that, in these asymptotic conditions, the equilibrium 
and shear viscous pressure decay much faster than mc 2 n, the rest mass 
energy. Therefore, the solutions behave asymptotically like dust solutions 
with the same sign for the function K. 

9.3 FLRW subcases. 

If together with spherical symmetry, we assume homogeneity by imposing 
"homogeneous" initial conditions: no = no and To = To, where no, To are 
positive constants, we have: 

M= mc2floY \ t7o = ^4 (48a) 
3 3(7-1) 

If we also assume K in (47b) to take the form: K = feo^o j we can re-write 
(21) as 



20 



. 2 K 

y = - 

y 



n + 



y 



3(7-1) 



, 2 - M _ Uo Y 
-koC, M=y3, W =y3> V= % ( ' 



so that its solutions will be of the form y = 
separable function given by Y — (R/ Ro)Yo, 
R(t ). From (47b) we have: 1/h 2 = 1 - fe F 2 
metric 



fit — to), and Y = $ is a 
where i? = i?(t) and Ro = 
, and (46) becomes a FLRW 



as 2 = -c 2 dt 2 + (Jp) [dr 2 + Y ( 2 {r) (d6 2 + sin 2 9d<j> 2 )] (48c) 

where Y ' — yjl — k~oY£, so that: Yo = (r, sin r, sinh r) for feo = (0, 1, — 1), 
Inserting Y = (i?/-Ro)Y into (22)-(29), we obtain cr = P = 0, 6 = 
3R lt /R, together with 

_ fRo\ 3 fR \ 3 [.->,« fRo\ 3 ^~ iy 

/ T? \ 3(7 — 1) / D \ 37 

fcT = fcT (^J , p=( 7 -l)f7 (^) (48d) 

and so, the source of this FLRW spacetime is a perfect fluid (P = 0) 
satisfying the generic equation of state (4) ' 26 ' . This particular case is the 
equilibrium limit of the solutions, characterized by fluid models of ideal 
gases following quasi-statical reversible processes along all fluid worldlines. 
Although the dust subcase can also be considered an equilibrium limit, it 
is a singular type of equilibrium since T = 0. 



10 Matching along a spherical boundary. 



We illustrate in this section the fact that the solutions can be smoothly 
matched along a spherical interface, either to a spherical or FLRW sub- 
cases, or to the Vaidya spacetime' 11 ' 12 ' . Since spherical symmetry follows 
by specializing the functions A,Bi,B2,C to the values given by (47), and 
since these functions determine v, no, To and are entirely arbitrary, it is 
possible to prescribe them so that spherical symmetry is smoothly reached 
for a given arbitrary r = r B . This requires the following boundary condi- 
tions associated with (47) on A,Bi,Bi,C 



Si(r B ) = S 2 (r B ) = 0, B' 1 ir B ) = B' 2 ir B ) = 
4A(r B ) = C(r B ) = const. A'(r B ) = C(r B ) = 



leading to 



v(r B ,0,4>) = l, u'(r B ,e,<t>)=0 
Y B =Y B (t) = Y{t,r B ,0,4>) = *{t,r B ), 



0) = 







V" 




V 


Y 


B 





(49a) 
(49b) 

(50a) 
(50b) 
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' dv~ 




~dv~ 


.de. 


r B 





= 0, 







' 8Y~ 


[de. 


r B 


d<j> 







(50c) 



where \\ r means evaluation along r — r B . Therefore, if conditions 
(49) hold, equations (50) guarantee that the hypersurface marked by 
(t, r B , e, <f>) is a world tube corresponding to the proper time evolution 
(t is proper time) of a class of 2-spheres labelled by r — r B , with proper 
radius Y B . Notice that this hypersurface has spherical geometry even if 
the geometry of the regions < r < r B and r > r B is not spherically 
symmetric. 

Assuming (49) and (50) to hold, we can consider the hypersurface 
r — r B as a spherical "boundary" separating two regions: (I) and (II) 
respectively marked by the ranges: < r < r B and r > r B . A situation 
of physical interest follows by looking for an "exterior" solution for (46), 
which means considering the system formed by: (I), the general solution 
described by (46), the "interior", and: for (II), a suitable spherically sym- 
metric spacetime, the "exterior" . Since each region (I) and (II) is a chunk 
of a different spacetime, we have actualy the problem of a matching be- 
tween two spacetimes. The conditions for smoothness of such matchings 
are the continuity along the matching boundary of the metric and the ex- 
trinsic curvature (K a b = —n a - c hi, where n a is a unit vector normal to the 



matching surface). If the same coordinates (t, r, e, 
(I) and (II), these matching conditions are ' 27 ' 



are used for regions 



9 t j 



(ii) 
9ij 



0. 



Ki 



9ij,r 



(51) 



where — (t,e,(f>) and [] B means evaluation at r = r B . It is straight- 
forward to prove that if the exterior spacetime is a particular case of (46) 
then (49)-(50) are equivalent to (51). We shall consider below the cases 
when region (II) is either one of: (a) a FLRW subcase, (b) a spherically 
symmetric but non-FLRW subcase, and (c) Vaidya spacetime. 



10.1 Matching with a FLRW or a spherical sub- 
case. 

If region (I) is characterized by a general case of (46) where A, Bi, B2, C 
comply with boundary conditions (49), then (50) and (51) also hold. If for 
the whole of region (II) we have: n = n (r B ) = n and T = T (r B ) = T , 
where no, To are positive constants, then region (I) is smoothly matched 
to a FLRW cosmology. However, the latter is not an arbitrary FRLW 
cosmology, but the particular case characterized by (48) . It is possible to 
reverse the roles so that a chunk of a FLRW spacetime and the general 
solution (46) respectively occupy regions (I) and (II). Matchings of this 
type, especially those in which the FLRW spacetime is the exterior, are 
theoretically interesting as models of inhomogeneities within a homoge- 
neous background. The possibility of performing such matchings has been 
reported' 25 ' for the perfect fluid Szekeres-Szafron solution (re-interpreted 
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as a mixture). The result presented here is a direct generalization to a 
viscous source. 

It is also possible to match a general solution (46) with spherically 
symmetric (but non-FLRW) subcase, occupying either region (I) or (II). 
Again, boundary conditions (49)-(51) must be satisfied and (47) must 
hold in the region occupied by the spherical subcase of (46). However, 
now T = To(r) and no = no(r) in that region. 



10.2 Matching with a Vaidya exterior. 

We shall examine the matching, along r — r B , between a general solution 
(the interior) and the Vaidya solution (the exterior). The Vaidya space- 
time is known' 11 ' 12 ' to be the exterior solution for spherically symmetric 
imperfect fluid sources and it is regarded as a radiating generalization of 
Schwarzschild, since its source is a null dust (or pure radiation) type of 
electromagnetic field. It is usually described by the metric 



ds 2 = — 



i 2m(v) 



Y 



c 2 dv 2 - 2c dv dY + Y 2 [d9 2 + sin 2 9d(j> 2 ] (52) 



where the null coordinate v is a "retarded time" . Let the interior be a 
general solution characterized by (46) with A,Bi,Bi,C and no, To satis- 
fying the boundary conditions (49)-(50) or (51) at the boundary r = r B , 
so that this boundary is a spherical region and the exterior r > r B is 
characterized by 

(Y'Y 

ds 2 = -c 2 dt 2 + - K _ ^ dr 2 + Y 2 [d6 2 + sin 2 9dcj> 2 ] (53) 

where Y = Y(t,r) (or, equivalently: Y = <f>) and K is given by (47b). If 
Y in (53) satisfies (21) with 



M = M B = r n Y 2 Y^dr, U = Uo B = P I^P^L, y = Yo B 

J r=0 J r — ^ 

(54a) 

then this exterior is the specific case of Vaidya's spacetime (in the coor- 
dinates t,r,6,4>) with 



2rh(v) = k 



M B +U 0[ 



V \ 3(7-1) 



B I y 



(546) 



where Y B = Y(t,r B ), Yo B = Yo(r B ), and the relation between v, t and 
Y B at the boundary is given by 



±1 



dv 



dY 



^l-K B ci±Y B 
±1 



b Y B [y/l-K B C*±Y B 



(55a) 
(556) 
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where Y B is (21) with Y = Y B and M, Uo,Yo given by (54a), or equiv- 
alently: Y B — 2m/Y B — K B c 2 . The integration of (55b) provides the 
trajectory of the boundary in the (v,Y) plane of Vaidya spacetime (52). 
For r > r B , the transformation relating the coordinates (t, r) in (53) with 
(v,Y) in (52) are 

Ydt + Y'dr = dY (56a) 

= cdv (566) 



VI - Kc 2 ± Y VI - Kc 2 [VI - Kc 2 ± V] 

where Y is (21) with M, Uo,Yo given by (54a) but K = if(r) arbitrary. 
By evaluating the corresponding mixed derivatives, it is possible to verify 
that (56) are well defined coordinate transformations. 

An interesting result is the fact that the the matching between a gen- 
eral solution and Vaidya spacetime is possible, without assuming that 
the interior region is spherically symmetric. This occurs also' 28 ' in the 
matching between Szekeres dust solution and Schwarzschild (the partic- 
ular case Uo = of the matching examined here). Another interesting 
point follows from (54b): the time dependent "Schwarzschild mass" con- 
tains the contribution from the rest mass (M) and internal energy (U) of 
the spherical region inside the boundary. If this boundary is expanding, so 
that: Y B /Yo g ^> 1, this mass function decreases tending asymptotically 
to a constant Schwarzschild mass, and indicating that the internal energy 
contribution to the total mass-energy of the spherical interior is radiated 
away until only the rest mass contribution remains. 



11 Concluding remarks. 

We have found a new class of exact solutions of Einstein's field equations, 
characterized by the Szekeres-Szafron metric of class I, whose source is a 
fluid with shear viscosity. The solutions admit, in general, no isometries 
and generalize well known solutions with dust and perfect fluid sources. 
However, their main appeal stems from the fact that the equilibrium state 
variables (equations (22)-(24)) satisfy a physically meaningful equation of 
state: that of a relativistic ideal gas in the following important limits: (1) 
ultra-relativistic, (2) non-relativistic and (3) a binary "matter-radiation" 
mixture of gases of the type (1) and (2). Since the process of integration 
of the field equations does not involve the shear viscous stress (equations 
(11) and (25)), we have provided the conditions constraining this quan- 
tity, and relating it to equilibrium variables and to an entropy production 
law, in terms of of the equations of Causal Irreversible Thermodynamics 
with phenomenological coefficients related to those obtained from Kinetic 
Theory for the Maxwell-Boltzmann gas and the "radiative gas" . Together 
with the conditions for thermodynamic consistency (equations (41)), we 
have given a regularity condition that prevents the emergence of unphysi- 
cal singularities (equation (42)), as well as examined the fulfilment of the 
weak, dominant and strong energy conditions (equations (43)-(45)). 
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The solutions become fully determined once the spatial gradients Y' /Y 
are explicitely computed from the integrals (31)-(33) of the quadrature 
(30) that follows from the Friedmann-like evolution equation (21). The 
temperature and particle number density, no and To, at an initial hyper- 
surface marked by an arbitrary t = to, together with a function related to 
the extrinsic curvature of this hypersurface (K), emerge naturally as ini- 
tial value functions. Once a specific solution is determined for any given 
form of these functions, it is necessary to test the fulfilment of regular- 
ity and energy conditions, as well as the conditions for thermodynamic 
consistency. This task is beyond the scope of the present paper and is 
already being undertaken in separate papers, dealing with the classical 
ideal gas' 10 ' and the matter-radiation mixture ' 23 ' . 

We have also looked at the matching conditions with FLRW and 
Vaidya spacetimes along a spherical boundary. Such a boundary can 
always be defined, even if the enclosed "interior" region is not spherically 
symmetric. The case of matching to a Vaidya exterior metric deserves es- 
pecial attention, since this solution is the usual "exterior" metric of local- 
ized (spherical) objects whose sources are imperfect fluids, and there is a 
large body of literature concerned with these models. These solutions are, 
either shear-free fluids with heat conduction' 29 ' , or fluids with anisotropic 
pressures' 30 ' (but without associating this anisotropy with viscosity). In 
many cases, especially those without shear, these imperfect fluids do not 
obey physically meaningful equations of state. The imperfect fluid sources 
in this paper, having a consistent thermodynamical interpretation, allow 
one to examine the geometry and dynamics of a more realistic source ex- 
panding or collapsing in a Vaidya spacetime. The introduction of thermal 
phenomena and radiative processses into a collapsing/expanding localized 
source might lead to interesting generalizations of the so called "Swiss 
Cheese" models (dust sources matched to vacuoles with Schwarzschild 
geometry, see chapter 3 of reference [11]). 

Those specific solutions complying with regularity and energy con- 
ditions and with thermodynamic consistency have an enormous poten- 
tial as models in applications of astrophysical and cosmological interest. 
Consider, for example, the following possibilities: (1) Cosmological voids. 
There is a large body of literature using LTB (Lemaitre-Tolman-Bondi) 
dust solutions to model the evolution of great voids (see chapter 3 of [11] 
for a review). Since the solutions in this paper generalize LTB dust solu- 
tions, allowing for a consistent decription of thermal phenomena, it should 
be possible to generalize previous work in order to be able to consider 
those cases in which the presence of thermal motions (i.e. associated to 
some kind of dark matter) yields a non-vanishing pressure. (2) Structure 
formation in the acoustic phase. Again, there is a large body of literature 
on the study of acoustic perturbations (see for example reference [13]) in 
relation to the Jeans mass of surviving cosmological condensations. Equa- 
tions of state identical to (5b) are often suggested in this constext' 3 '"' 4 ''' 16 '. 
Since practically all work on this topic has been carried on with perturba- 
tions on a FLRW background, the exact solutions derived and presented 
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here, especially the matter-radiation mixture, are idealy suited as an al- 
ternative treatment for this problem. (3) Effect of inhomogeneities on the 
microwave cosmic background. Recent papers' 31 ' have applied numerical 
techniques (in the non-linear regime) to explore possible anisotropies in 
the microwave background due to photons crossing inhomogeneities, the 
latter being modelled by a LTB dust solution. It is worthwhile considering 
the extension of this work to encompass a solution whose source can be a 
thermodynamically consistent photon gas or mixture of a photon gas and 
a non-relativistic ideal gas. These and other applications are worth to be 
undertaken in future research efforts. 
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